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Using Langevin simulations we examine driven colloids interacting with quenched disorder. For 
weak substrates the colloids form an ordered state and depin elastically. For increasing substrate 
strength we find a sharp crossover to inhomogeneous depinning and a substantial increase in the 
depinning force, analogous to the peak effect in superconductors. The velocity versus driving force 
curve shows criticality at depinning, with a change in scaling exponent occuring at the order to 
disorder crossover. Upon application of a sudden pulse of driving force, pronounced transients 
appear in the disordered regime which are due to the formation of long-lived colloidal flow channels. 
PACS numbers: 82.70.Dd,64.60.Ht 



Colloidal crystals are an ideal system in which to study 
the general problem of ordering and dynamics in 2D 
10-^ , since the particle size permits direct imaging of the 
particle locations and motion. A considerable amount of 
work has been conducted on the melting of 2D colloidal 
crystals in the absence of a substrate ■ In addition, 
a number of experimental and theoretical studies have 
considered colloidal crystallization and melting in 2D sys- 
tems with periodic ID |^,|| and 2D substrates where 
a rich variety of crystalline states can be stabilized. 

Colloid crystals are also ideal for studing the ordering 
and dynamics of an elastic media interacting with ran- 
dom substrates, a problem that is relevant to a wide vari- 
ety of systems such as superconducting vortices, Wigner 
crystals, and charge-density waves (CDWs). Open issues 
include the nature of the dynamical response to applied 
forces, as well as whether an order to disorder transition 
occurs as the strength of the random substrate increases. 
Recently, LeDoussal and Carpentier have theoretically 
investigated the effects of quenched disorder on the or- 
der and melting of 2D lattices and find a sharp crossover 
from the ordered Bragg glass (where defects are absent) 
to a disordered or molten state |^ . They predict that the 
depinning threshold increases at this crossover due to the 
softening of the lattice, which allows the particles to bet- 
ter adjust to the substrate. A similar mechanism could 
account for the peak effect observed in vortex matter in 
superconductors [p|-p^, in which the depinning thresh- 
old rises dramatically when the applied magnetic field is 
increased. In low temperature superconductors, where 
the vortices are fairly stiff so that their behavior can be 
considered as effectively 2D, recent small angle neutron 
scattering experiments have shown that the peak effect is 
associated with a sharp disordering or melting transition 

In addition to static properties, the dynamics of elastic 
media interacting with quenched disorder in 2D is a topic 
of intense study. In the disordered region the driven sys- 
tem may break up into pinned and flowing regions, as ob- 
served in experiments and simulations ]l6|-p^ of vor- 
tices in superconductors. Conversely, for weak substrate 



disorder the elastic media is defect free and undergoes 
elastic depinning, in which the particles keep the same 
neighbors as they move. Fisher predicted that elastic de- 
pinning would show criticality |p3|| , and that the velocity 
vs. force curves would scale as v = {f — fc)^ , where /c 
is the depinning threshold force. This scaling has been 
studied extensively in 2D CDW systems where /? = 2/3 
p^-p6|. It is, however, not known whether this exponent 
occurs in other systems undergoing elastic flow. Another 
intriguing dynamical phenomena is the pronounced tran- 
sient behavior exhibited by vortices under a sudden ap- 
plied current pulse at magnetic fields near the peak ef- 
fect regime 12 . Due to surface barrier effects in the 
vortex sample, it is not clear whether these transient ef- 
fects relate to the plasticity of the vortex dynamics or to 
contamination of the vortex lattice by disorder from the 
sample edges [Q . Recently, Pertsinidis and Ling Q have 
studied colloids in 2D driven by an electric field and inter- 
acting with a disordered substrate. They observe plastic 
depinning with filamentary or river-like flow of colloids 
and a velocity- force curve scaling with /3 = 2.2. Under a 
pulsed drive the system shows very long time transients 
that fit to a stretched exponential. 

Motivated by the recent colloidal experiments as well 
as the pulse drive experiments in vortex matter, we have 
conducted Langevin simulations of colloidal particles in- 
teracting via a Yukawa potential in 2D systems with ran- 
dom disorder. In simulation, the strength of the disorder 
can be carefully tuned, which is difficult to achieve in ex- 
periments. In addition, the initial conditions of the col- 
loidal arrangements are easily controlled, whereas in ex- 
periments, defects generated in the colloidal lattice dur- 
ing preparation may become frozen in by the disorder. 

Our simulations show that for weak substrates the col- 
loids form an ordered triangular array which depins elas- 
tically without the generation of defects. For increased 
substrate strength, there is a sharp crossover to a dis- 
ordered phase where the colloids depin plastically into 
riverlike structures. This crossover is accompanied by a 
sharp increase in the depinning threshold, analogous to 
the peak effect phenomena in superconductors. We 
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tex simulations produce results similar to the parabolic 
pinning |la-n8| . The equation of motion for colloid i is 



(d) 



FIG. 1. (a) Depinning force fc vs pinning strength fp. 
Inset: Corresponding derivative dfc/dfp. (b) Percentage of 
non-sixfold coordinated colloids, Pd- (c) Delaunay triangu- 
lation of colloid positions at depinning for fp = 0.12; (d) 
fp = 0.25. Filled circles indicate defect sites of particles which 
are not sixfold coordinated. 

find scaling of the velocity vs. applied drive with an 
exponent of /3 = 0.67 in the elastic regime, in agreement 
with studies in 2D CDW's. In the plastic regime we find 
/3 = 1.97, close to the value observed in experiments [Q. 
In the disordered region, long time transients that fit to 
a stretched exponential occur in response to a sudden 
applied drive pulse, as also observed in the colloidal and 
vortex experiments. 

The colloids are simulated using Langevin dynam- 
ics in 2D The colloids interact via a Yukawa 
or screened Coulomb interaction potential V{rij) — 



il)exp(- 



Here Q is the charge of 



the particles, 1/k is the screening length, and is the 
position of particle The length of the system is mea- 
sured in units of the lattice constant ao and a screening 
length oi K = 2/ao is used. The quenched disorder on the 
substrate is modeled as randomly placed parabolic traps 
with radius rp < qq and a maximum force fp. This same 
type of model for pinning has been used previously to 
model quenched disorder in superconducting vortex sys- 
tems fly] . Other types of pinning potentials used in vor- 
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Here = Viy(rij) is the interaction force from 

the other colloids, fp is the pinning force, is a randomly 
fluctuating force due to thermal kicks, and is the force 
due to an applied driving field. We start the system at a 
temperature above melting as determined from the dif- 
fusion, and gradually cool below to T/Tm = 0.4 for 
most of the data presented here. To measure velocity 
V vs force curves, care must taken to average over sub- 
stantial amounts of time to avoid transient effects. The 
driving force is increased from zero by small increments 
and the velocity is averaged for 5 x lO'' time steps at 
each increment, with typical simulations running for 10^ 
time steps. In this model we do not take into account 
hydrodynamic effects or long-range attractions between 
colloids. 

In Fig. 1(a) we show the depinning force fc vs. sub- 
strate strength fp from a series of simulations. For 
fp < 0.18 the depinning force increases as a power law, 
fc on /p~^'^*°'^. To compare the depinning force to the 
order in the system, in Fig. 1(b) we show the the per- 
centage of defects or non-six fold coordinated particles 
Pd as calculated from a Delaunay triangulation. This 
measure indicates that the colloidal crystal is in an or- 
dered state {Pd = 0.0) for fp < 0.18, and that there is a 
crossover to a disordered state {Pd ^ 0) at fp = 0.18. In 
Fig. 1(c) we show a representative Delaunay triangula- 
tion for the ordered state where there are no defects but 
small distortions in the particle positions can be seen, 
and in Fig 1(d) we show the disordered state where de- 
fects are present. The crossover to the disordered state 
coincides with a rapid increase in the depinning force as 
seen in Fig. 1(a) and in the inset of Fig. 1(a) which shows 
a peak in dfc/dfp at the crossover. This behavior is con- 
sistent with the recent experiments in superconductors 
which find that at the peak effect there is an increase in 
the pinning with a simultaneous disordering of the lat- 
tice |l^. For fp > 0.2, the depinning scales as fc oc /p, 
as expected for the single particle pinning regime. The 
sudden increase in the depinning force results from the 
fact that the defected colloid lattice is much softer than 
the ordered lattice, allowing the colloids to adjust their 
positions to accommodate to the optimal pinning sites. 
We have also investigated this transition for different col- 
loidal densities and disorder strengths. For increasing T 
the order to disorder transition is shifted to lower val- 
ues of fp. We have also investigated finite size effects for 
increasing system sizes, and find that the order-disorder 
crossover shifts only a small amount before saturating, 
while the sharpness of the transition persists with in- 
creased system size. 
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FIG. 2. Colloid positions (black dots) and trajectories 
(lines) for (a) plastic flow regime {fp — 0.25), and (b) elastic 
flow regime {fp — 0.12). 

It is beyond the scope of this paper to determine 
whether the order to disorder crossover is a first order 
transition. Although the sharpness suggests a possible 
first order transition, Carpentier and LeDoussal show 
that for 2D systems with quenched disorder, a sharp dis- 
ordering crossover, rather than transition, occurs [||. In 
addition, a first order transition is not expected since the 
Bragg glass in 2D has been shown to have dislocations 
on large scales at all temperatures. The distance between 
these dislocations can be arbitrarily large |28|. 

In Fig. 2 we show that the order-disorder crossover co- 
incides with the onset of plastic flow above depinning. In 
Fig. 2(b) the elastic colloid flow is shown for fp = 0.12 
above depinning {fd/fc — l-l)- Here each colloid keeps 
the same neighbors as it moves. In Fig. 2(a) the inhomo- 
geneous or plastic colloidal flow is shown for fd/fc = 1-1 
for fp — 0.25. Here, only a portion of the colloids are 
moving at any one time, and the motion occurs in chan- 
nels or rivers between pinned regions. The colloid ve- 
locities show a bimodal distribution in this regime, split 
between the stationary and moving colloids. In addition, 
the channels seen in Fig. 2(a) are not static but change 
over time, so that any one colloid is only temporarily 
trapped in a pinning site. These features of the plastic 
flow are in agreement with observations in colloidal ex- 
periments [0 and in vortex simulations of the strongly 
pinned regime p6|-[l9|. 

In order to correlate the different types of flow observed 
in Fig. 2 with properties of bulk measurements, we show 
in Fig. 3 the scaling of the velocity vs driving force. For 
elastic depinning in the ordered regime [Fig. 3(a)], the 
velocity vs force curves fit well to v = {fd — fcY with 
/? = 0.66 ±0.02, as illustrated in Fig. 3(b). These results 
are in good agreement with theoretical predictions p^ ] 
and simulation results |25 2^ for CDW's in 2D where 
the depinning is elastic. In contrast, in driven 2D vortex 
matter, Higgins and Bhattacharaya ||l0|] found an expo- 
nent of /3 = 1.2 below the peak effect where elastic flow 
is expected to occur. This may be due to the effects of 
surface barriers disordering the lattice. In Fig. 3(c,d) 
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FIG. 3. Velocity v vs applied drive fd for (a) elastic regime 
fp — 0.08. (b) Log-log plot of V vs {fd — fc) from (a); line 
indicates fit to /3 = 2/3. (c) v vs fd for plastic depinning 
fp = 0.25. (d) Log- log plot of V vs {fd — fc) from (c); line 
indicates fit to /3 = 2.2 

the V — fd scaling for the plastic regime is presented, 
with f3 = 1.94 ± 0.03, close to the value of 2.2 found in 
the colloid experiments 0. For larger system sizes, we 
find that the scaling region is expanded but the exponent 
is unchanged. Other studies in the plastic flow regime 
found /3 = 5/3 for electron flow in metallic dots ||2^ and 
(3 = 2.22 for vortex flow in Josephson-junction arrays 

ill- 

In Fig. 4 we show the response of colloids prepared 
in an ordered state to the application of a sudden pulse 
of driving force of different strengths in the plastic flow 
regime. Since the pulse strength is chosen to be below the 
depinning threshold value fc, the initial colloid velocity is 
high, and then gradually decreases. We find that a sim- 
ple functional form cannot be fit to the curves. Instead, 
we use a stretched exponential fit as performed in exper- 
iments 0: v{t) = voexp{-{t/to)°')+vi. The values of to 
and a depend on the magnitude of the drive. For the pa- 
rameters investigated here, a falls between 0.08 and 0.4, 
in agreement with the values found in experiment. A sim- 
ilar stretched exponential decay was also found in vortex 
matter for the transient response to pulses |l2| . We find 
that in the long time limit, the colloid flow occurs only 
through a few long-lived channels. In the elastic regime, 
the decay of v is much faster and fits to an initial pure 
exponential with the velocities going to zero. We note 
that in the elastic regime the colloids move less than a 
lattice constant after a pulse is applied, whereas in the 
plastic regime, colloids in the moving channels can move 
the entire length of the system. For increased system 
sizes, the transient times are enhanced in the plastic flow 
regime but are unchanged in the elastic regime. The long 
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FIG. 4. Average transient velocity per colloid V/Nc vs 
time in the disordered region for /p = 0.25 at a sudden applied 
drive of (from top to bottom) ,fd/ fc ~ 1.2, 1.0, 0.8, 0.6, and 
0.3. A stretched exponential A exp(— (t/to)") + vi can be fit 
to all the curves. Bottom curve: a stretched exponential fit for 
the fd/fc = 0.3 case, with V/N^ = 0.0089 exp(-(t/400)^) +0. 

lived transients in the plastic regime are responsible for 
the very slow velocity-force sweep necessary to measure 
an accurate depinning threshold. This sweep-rate depen- 
dence is also consistent with the experimentally observed 
sweep rate dependent critical currents in the peak regime 
where slow rates produce larger measured critical 
currents. 

To summarize, we investigated the behavior of 2D col- 
loids interacting with random disorder using Langevin 
simulations. For weak disorder the colloids form an or- 
dered lattice which depins elastically and shows critical 
scaling in the velocity vs force curves, with /3 = 0.67, in 
agreement with studies of 2D CDW's. For increasing dis- 
order strength, we find a sharp crossover to a disordered 
state, accompanied by a sharp increase in the depinning 
force, analogous to the peak effect observed for vortex 
matter in superconductors. In the disordered region, the 
colloids depin inhomogeneously into fluctuating channels 
and the v — f curve scaling gives /3 — 1.97, in agreement 
with experiments. In the disordered How regime, pro- 
nounced transients occur in response to a sudden pulse, 
with the late time dynamics determined by a few long 
lived channels. 
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